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Direct-imaging searches for planets reveal wide orbit planets amenable to spectroscopy, and their atmospheres represent an important 
comparison to the irradiated atmospheres of Hot Jupiters. Using AO integral field spectroscopy of 2M 1207 b, the shape of the con- 
tinuum emission over the J, H, and K bands from the atmosphere of this young, planetary mass companion is measured in order to 
compare with atmospheric and evolutionary models, and objects of similar temperature in young clusters and the field. The 2M1207 b 
spectrum has the highest spectral resolution (R~300-1500) and largest wavelength coverage, including the first /-band spectrum, 
for this benchmark object. The high signal-to-noise of the data allow a clear identification of signatures of low surface gravity, and 
comparison with a grid of AMES-Dusty models reveals a best-fit effective temperature of T sff = 1 600 K with a preferred surface gravity 
of log g=4.5 . The ./-band flux is depressed relative to nearly all L-type objects, and the detailed shape of the absorption features across 
the //-band exhibit differences from the model predictions. The possible origins of 2M1207b and its low luminosity are examined 
with the new data and analysis which suggest that extinction from a disk with large grains is a viable scenario and is preferred over 
scatttering off an optically thick disk. The 2M1207 b spectrum presents an important comparison for the types of features which may 
be present in upcoming spectra of the atmospheres of planets imaged in orbit around stellar primaries. 
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1. Introduction 

Over 400 extrasolar planets are currently knowrQ with most sys- 
tems detected indirectly in orbits closer than that of Jupiter. For 
indirectly detected planets, the rare combination of a favourable 
geometry resulting in a transit and a very bright host star for 
high signal-to-noise data is required to obtain a planetary spec- 
trum either in transmissi on (|Charbonneau et al.|2002"| ) or emis- 
sion ( |Charbonneau et al.||2T}05||Deming et al.|2005^ Although 
only two transiting planets have been investigated spectroscopi- 
cally, transmission studies have shown the possibility of methane 



gation of the early evolution of planetary atmospheres. The first 
spectrum of HR 8799 c in the L-band suggests differences in the 
continuum shape from theoretical models (Jan son et al.||2010 



absorption ( |Swain et al.|2008)> or Rayleigh scattering (Sing et al. 
2008) and a haze layer ( Pont et al.|2008|l, and emi ssion studies 
have observed the presence ^Grillmair et al.|2008) > and absence 
( |Rhee et al.|2007] > of water band absorption in different systems. 

The atmospheres of directly-imaged planets represent a com- 
plementary case study to the transiting planet spectroscopy, un- 
contaminated by the physical effects of extreme proximity to the 
host star. In some cases, spectroscopy of directly imaged planets 
is feasible. The first confirmed directly-imaged planets orbiting 
normal stars ( |Kalas et al.|2008t|Marois et aL|2008|) and the first 
imaged planetary mass companion (Chauvin et al. 2004) present 
the opportunity for spectroscopic analysis of young planet at- 
mospheres. For technical reasons, targets observed in imaging 
planet searches are typically younger than stars searched for 
planets indirectly, and this selection criterion enables an investi- 



possibly due to non-equilibrium CO/CH4 chemistry and indicat- 
ing the importance of measuring the atmosphere over an even 
larger wavelength range. This paper presents the results of the 
highest resolution J, H, K spectrum of 2M1207 b which should 
provide a high signal-to-noise comparison spectrum to planets 
imaged in orbit around stellar primaries. 



2. Target 

The target for this study is one of the lowest mass compan- 
ions imaged to date ( [Chauvin et~aT1[2004l |2005] ). 2M 1207 b 
is confirmed as a companion based on multi-epoch imaging 
( |Song et al.|2006 l and, considering the distance of 52.4+1.1 pc 
( |Ducourant et al. 2008 ), orbits the host brown dwarf at a pro- 
jected separation of 46 AU, somewhat wider than the Sun- 
Neptune distance. The companion has apparent magnitudes 

~ |2007) , #=18.09+0.21, and 



of 7=20.0+0.2 (Mohanty et al 



/T s =16.93±0.11 ( |Chauvin et al.||2004|). Combining the system 
age and theoretical cooling tracks ( Chabrie r et al.|200Tj{ [B araffe 
|et al.|[2002] |2003) > with the target photometry ( |Chauvin et al. 



* Based on observations obtained at the Paranal Observatory, Chile 
for ESO program 078.C-0800(B). 
1 www.exoplanet.eu 



120056 indicates a mass of ~5 M Jup , whereas with the spectro- 
scopically estimated effectve temperature, a mass of 6-10 Mj up 
is predicted (Moha nty et al.|20 07). The masses depend critically 
on the system age, estimated to be 8 Myr for 2M1207 from its 
membership in the TW Hydra association (Gizis 2002); the un- 
usually well-determined age of this low mass companion elimi- 
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Table 1. SINFONI observations 



Source 


Grating 


Scale 


NDIT 


DIT 


#exp. 


Standard (HIP) 


T cff (K) 


Dates 


Project ID 


2M1207A+b 


H+K 


100 


1 


300 s 


24 


54930, 62566 


29000, 15000 


28 Jan 07, 07 Feb 07 


078.C-0800(B) 




J 


100 


1 


600 s 


6 


59951 


16400 


22 Feb 07 


078.C-0800(B) 



nates some degeneracies in model comparisons. Given the large 
difference in mass ratio of the 2M1207 system compared to the 
HR 8799 and Fomalhaut systems, their formation mechanisms 
may be different, with 2M1207 representing the tail of the bi- 
nary distribution and the A-star planets examples that formed 
in a disk via core accretion (e.g., Pollack et al.||1996*] l or grav- 
itational instability (e.g., Boss|1997| >. Despite the possibility of 
different origins, 2M1207b is a benchmark object for the study 
of young planetary atmospheres. 

3. Observations 

The observations were taken with SINFONI, an AO-equipped 
integral field spectrograph mounted at the Cassegrain focus of 
the VLT ([Bonnet et al.|[2004l |Eisenhauer et~aLl[2003) |Thatte| 
|etaT|l998| > in program 078.C-0800(B) (PI: N. Thatte). Details 
of the observations are given in Table[TJ The primary object 
was used for the AO correction. Of the four gratings available 
in SINFONI, the / (R-2000) and H+K (R-1500) were used, 
which yields higher spectral resolution than previous results 
( |Chauvin et al.|2004| |Mohanty et al.|2007| l and the first /-band 
spectrum. The 100 mas spatial scale was employed as it is fine 
enough to separate the pair, but coarse enough to concentrate the 
flux of the faint companion. The field-of-view at this scale is 3'.'0 
x 3'.'0 and included both components. 

The observing sequence included a set of exposures on the 
target at several offset positions on the array. Observations of 
an early type star followed each target sequence to correct for 
telluric features and the instrument response. Additional calibra- 
tions were obtained to measure the wavelength scale, the detec- 
tor dark current, distortion, and quantum efficiency variations. 



4. Data Analysis 



The data were reduced with the Gasgano (Gebbinck et al. 
|2007|l implementa tion of the SINFONI data reduction pipeline 
( jModigliani et al.||2009| |Dumas et aT|2007| l. With a set of de- 
tector calibration files, corrections for the dark current, linearity 
response, and bad pixels were measured and applied to all data 
files. Two effects of the instrument optics were also corrected 
- the distortion of the pixel scale across the field and the non- 
linear dispersion of the gratings. After the detector and optics 
processing was completed, the array data were converted into 
data cubes consisting of a stack of images, each covering a slice 
of wavelength. Across the entire wavelength range, a stack of 
images were constructed, composed of ~700-2000 wavelength 
slices, depending on the band. Binning was required for the J- 
band data. 

The final steps applied to the data cube involved sky sub- 
traction, combining the offset positions, radial subtraction of the 
host object, and tracing the path of the companion spectrum 
through the data cube. The targets were observed at different 
positions on the array in an ABBA pattern sampling two areas 
on the array. The cubes with the target in one portion of the array 
were averaged and then used as the sky to subtract off the indi- 
vidual cubes at the opposite position. Next, the individual sky- 
subtracted cubes were combined, centroiding on the primary. 




Fig. 1. The J,H,K wavelength collapsed images before and after 
the alignment and subtraction process. For the H— and K- band 
portions of the data, the companion is very distinct; although 
the /-band image is lower signal-to-noise, the companion is 
still clearly evident, enabling the spectrum to be measured in 
the binned data cube. 



The individual wavelength planes were shifted and aligned on 
the primary such that the final path of the companion through the 
data cube was a straight trajectory. Due to the small angular sep- 
aration between the primary and secondary, the halo of the pri- 
mary was fit with a radial profile and subtracted from each wave- 
length plane. Images at /, H, and K, collapsed in wavelength 
over the entire band, before and after the subtraction process are 
shown in Fig.[TJ The measured FWHM of the target was stable 
at ~3.5 pixels. Following the radial subtraction, the companion 
spectrum was constructed by extracting the flux in an aperture 
of 4.5 pixel radius that encompasses the majority of the flux in 
a consistent manner. Given the stability of the primary FWHM 
throughout each wavelength band, variations of the Strehl ratio 
across a band seem minimal and should not systematically bias 
the extracted spectrum. 

Correction for telluric absorption and the detector response 
was accomplished by dividing the observed target spectra by 
those of early-type standard stars with interpolation across in- 
trinsic features. The shape of the standard was then removed 
by multiplication by a black-body at the appropriate effective 
temperature (see Table[T]i. The extracted / and HK spectra of 
2M 1207 b were flux calibrated using the / and /Ts-band pho 
tometry of |Mohanty et al.| ( |2007| l and |Chauvin et~aL] ( |2004 



respectively. The expected ratio of detected photons was ca 
culated by convolving the observed spectrum and that of Vega 
( |Mountain et al.|1985||Hayes|1985| > with the 2MASS responses 
curves (Carpenter 2001 1, assuming Vega to have zero magnitude 
in both filters. 

The flux calibrated near-IR spectrum of 2M1207 b is shown 
in Fig. [2] Over the /-band the signal-to-noise ratio is ~1 per res- 
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Fig. 2. The spectrum of 2M1207b covering the /, H, and K passbands (black lines) flux calibrated using the apparent magnitudes 
of Chauvin et al. ( 2004| l and Mohanty et al. (2007). Key features are indicated. 



olution element, increasing to ~13 across the //-band and ~24 
over the Zf-band. The depression of / due to atmospheric dust 
absorption, the two prominent water absorption bands between 
the near-IR passbands, and the CO absorption at the red end of 
the ZT-band are clearly evident. 

To estimate the effective temperature and surface gravity of 
2M1207b, we compared the flux calibrated JHK spectrum to 
the AMES -Dusty model grid ( |Allard et aLpOOT) with 7/ eff =500- 
4000 K and log g=3. 5-6.0. The model spectra were smoothed 
to the same resolution as the observations and interpolated to 
the same dispersion. Using the same least-squares statistic as 
Mohanty et al. (2007) and allowing the model flux to be scaled 
to find the best global fit, we identified the best fitting models. 
Calibrating the model spectra to the flux as would be observed 
from Earth requires a scale factor which depends on the radius of 
the object and its distance. Since we have flux calibrated the ob- 
served spectrum, we can extract an estimate of the radius implied 
for the closest match between each model and the observations. 



5. Results and Discussion 

The larger coverage in wavelength and higher signal-to-noise 
enables an investigation of the atmosphere and origin of this 
young, planetary mass companion. The flux calibrated spectrum 
of 2M1207 b is compared with previous observations of field and 
young low mass objects and with theoretical model atmospheres. 
The overall shape is significantly different to field L dwarfs, with 
a marked depression in the /-band as shown in Fig. [3] The H- 
band continuum exhibits a distinct triangular shape unlike the 
older, higher surface gravity objects, and the A"-band peaks at 
~2.3 [im (as suggested by Moha nty" et al.| [2007 ) rather than at 
~2.15 /mi as with field L dwarfs, and the region associated with 
the peak is more clearly defined in these higher resolution data. 

The spectrum of 2M1207 b is unlike field L-dwarfs spanning 
a range of spectral types, but the young age of the TW Hydra 
association (8 Myr) complicates the comparison. In Fig.|4] the 
spectrum of 2M 1207 b is plotted with several spectra of low 
mass members of Upper Sco ( |Lodieu et al. 2008| l, with spec- 
tral types ranging from M8-L2, and with the borderline brown 
dwarf/planet companion AB Picb (Bonnef oy et al.|2010 l. Even 
amongst young, low-mass objects in the 5 Myr Upper Sco re- 
gion, the 2M 1207 b spectrum appears distinct in the shape and 
relative fluxes of the continuum across the different bands. No 



Upper Sco member has a lower /-band relative to //-band and 
the triangular shape of the //-band is only beginning to develop 
by the L2 spectral type. The spectral sequence shown in Fig. [4] 
suggests that 2M 1207 b is one of the coolest young companions 
imaged. The most analogous object to 2M1207b is the compan- 
ion AB Pic b which also shows some of the key features present 
in 2M1207b. Since the /-band data required binning, the spec- 
tral resolution is not yet high enough to identify and measure 
absorption lines, though future observations with more signal- 
to-noise could investigate the depths of the Nal and KI lines. 

A series of spectral indices based on ratios of the flux of 
different regions of the continuum have been developed to char 
acterise the spectral type of cool objects ( |McLean et al.|[2003 



|Slesnick et al.|2004[ |Allers et al.|2007| >. We measured th ese for 
2M 1207 b and summarise the results in Table|2]for the lAllersl 



et al. ( 2007| > H2O and gravity sensitiv e indices, the|Sles nick et al. 



P004| l H 2 Q and FeH indices, and the |McLean et al.| ( |2003| l HjO 



indices. These suggest a spectral type in the broad range M8.5- 
L4, while the gravity sensitive index is consistent with that found 
for other young brown dwarfs. The discrepancy in inferred spec- 
tral types is not surprising given the unusual overall shape of 
the spectrum and highlights the importance of obtaining larger 
wavelength coverage of flux calibrated spectra in assigning spec- 
tral types and their uncertainties. 

The flux calibrated 2M1207 b spectrum with R-300 at / and 
R~1500 at HK is shown in Fig. [B ] along with the best fit to the 
model DUSTY spectrum ( jAllard et aT1|2001| l. Available mod- 
els with temperatures higher and lower by 100 K are also plot- 
ted in Fig. [5] The best fit, estimated by a least-squares fit to the 
grid of DUSTY models, and confirmed by visual examination 
is shown in the middle panel of Fig. [5] with / e ff=1600 K and 
logg=4.5. The uncertainty in temperature is clearly less than 
±100 K, but the ability to discriminate surface gravity variations 
with the broad morphology and flux level of the spectrum is less 
than for effective temperature. Consequently, logg=4.5 is only 
marginally favoured over log g=3.5 or 5.5 due to the uncertainty 
on the flux level of the /-band spectrum resulting from the pho- 
tometry uncertainty of 0.2™. While the location of the peak in the 
//-band is consistent with the models, the observed slopes on ei- 
ther side are not as steep as the model predicts, possibly due to 
difficulties modelling the H2 collisionally-induced absorption. If 
the precision of the /-band photometry could be improved, the 
model variations would allow surface gravity to be constrained 
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Fig. 3. The J,H,K spectrum of 2M 1207 b (black line) and sev- 
eral comparison spectra of field L dwarfs with spectral types LI 
(blue), L4.5 (green), and L8 (red). All spectra are normalised 
over the 1 .60-1 .65 fim range. 

Table 2. Spectral indices from Allers et al. (2007), Slesnick et 
al. (2004), and McLean et al. (2003). 



Index 


H 2 


grav. 


H 2 0-1 


H 2 0-2 


FeH 


Value 


1.17 


1.01 


0.660 


0.831 


0.877 


SpT 


L0 


low grav. 


L0 


L1.5 


M8.5 


Index 


H 2 0-A 


H 2 0-B 


H 2 0-C 


H 2 0-D 




Value 


0.593 


0.724 


0.733 


0.817 




SpT 


L2.5 


L3 


L0 


L4 





to ~0.5 dex, important to test one formation model. These phys- 
ical parameters are consistent with those measured previously 
( |Mohanty et"aiT|[2007 



though the inclusion of the /-band and 
the higher signal-to-noise of these data makes the fitting more 
secure. 

We did not compare the data to the COND models since 
they do not adequately reproduce the previous lower resolution 
H- and K-bwA spectrum ( |Mohanty et al.||20"07 >. In addition to 
the treatment of dust condensation and settling, the inclusion of 
non-equilibrium CO/CH4 chemistry can significantly alter the 
shape of the continuum emission (Fortney et al. 2008), particu- 
larly in the K- and L-bands. A suitable parameterisation of the 
eddy diffusion coefficient for such models may reproduce the 
shape of the observed spectrum, however, the flux level for the 
T e g = 1400 K (the hottest shown), K zz =4 non-equilibrium atmo- 
sphere shown in |Fortney et al. (2008 1 is still significantly higher 
than 2M1207b; for the range of K^, log§, and metallicity pre- 
sented for the Jeff = 1400 K case, the /T-band magnitude varies by 
up to 0.7 mag. 

As noted previously, to match the predicted spectral model 
fluxes to the observed photometry, either 2M1207b is underlu- 
minous or has a smaller than expected radius. The required ra- 
dius to explain the photometry is ~O.O52/?0, about one third 
of the radius of ~0. 16/?© from the DUSTY evolutionary model 
prediction for an object of / e ff=1600K at 5-10 Myr, and signif- 
icantly smaller than any predicted radius for a partially degener- 
ate object above 1 M Jup ( Chabrier et al.|2009 1. 

The origin of 2M1207 b and its underluminous flux is uncer- 
tain, and different scenarios have been proposed for the specific 
case of 2M1207b: an edge-on disk (Mohant y et al.|2007 1 and the 
result of a collision of protoplanets (Mamaj ek & Meyer||2007| ). 



\B Pic B 



2M 1207b 



M8 
L0 



1.2 1.4 1.6 1. 



2.2 2.4 



Wavelength ([im) 



Fig. 4. The J,H,K spectrum of 2M 1207 b (bottom) and several 
comparison spectra of young (5-30Myr) objects in Upper Sco 
of spectral types (from the top): M8, L0, L2, and -L0-L1 (AB 
Pic b). The spectra have been flux calibrated as in Figure 1, and 
then offset by an integer for clarity. 



More generally, the impact on the luminosity due to initial condi- 
tions and the early accretion history of young objects with ages 



comparable to 2M1207b, has been investigated (Marley et al 
2007 Baraffe et al.|2009 1. The protoplanet collision hypothesis 



predicts that the object is much smaller and should have a surface 
gravity of 3.0 rather than 4.5. The upper panel of Fig. [5] shows 
the data with DUSTY spectral models with three different effec- 
tive temperatures for a lower surface gravity of log §=3.5 (since 
log §=3.0 is not available). The data favour the higher surface 
gravity log §=4.5 value, though increasing the /-band flux by 
one sigma would make the data consistent with the log §=3.5 
model. A /-band magnitude difference with smaller uncertain- 
ties would be required to entirely rule out the collision scenario. 

For an initial investigation into the edge-on disk possibility, 
we compared the photometry and spectra with synthetic spec- 
tra ge nerated by the TORUS and MCFOST codes ( |Pinte et al. 
2009). A central object simulated by the r eff =1600 K DUSTY 
atmosphere was encircled with a disk 10 AU in radius, consis- 
tent with the physical separation of the 2M1207 pair. The disk 
was constructed with canonical values for flaring, surface den- 
sity power law, dust size distribution and chemistry. The inclina- 
tion was allowed to vary from edge-on to the ~60 degree limit 
suggested for the geometry of the primary disk. Within the lim- 
ited range of disk parameters searched, scattering alone could 
not explain the flux level of the central object since the predicted 
scattered flux was only a few per-cent, rather than the required 
~10%. For comparison, the T Tauri stars HH 30 and HK Tau 
with imaged edge-on disks (Stapelfel dt et al.|[T9 98) are fainter 
by 3-4 magnitudes relative to the photospheric levels, consistent 
with a scattered flux level of a few per-cent or less. Thus, among 
disk scenarios, grey extinction from a disk composed of larger 
grains than considered in the simulations is favoured over scat- 
tering off the surface of an optically thick disk; grain sizes larger 



than 4fim have been suggested as a possibility (Mohanty et al 
2007 ). Additionally, extinction of a hotter object is not able to 
explain the wavelength dependence of the photometry over the 
entire 0.9-3.6 fim range. 

In summary, the SINFONI spectrum of 2M1207 b represents 
the highest resolution, largest wavelength coverage spectrum of 
a young, planetary mass companion. The overall shape is dis- 
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Fig. 5. The spectrum of 2M1207 b covering the J, H, and K passbands (black lines). Middle panel: the best-fit DUSTY model with 
r e ff=1600 K and logg=4.5 in red along with the 1500K (green line) and 1700K (blue line) logg=4.5 models. Top and bottom 
panels: DUSTY models of the same effective temperature as the middle panel, but with surface gravity of logg=3.5 and logg=5.5, 
respectively. The observed spectrum has been flux calibrated using the apparent magnitudes of Chauvin et al. (2004) and Mohanty 
et al.| ( 2007] l. While the DUSTY models reproduce the general morphology, the flux scaling necessary implies a radius of ~0.052/?q 



at a distance of 53 pc, significantly smaller than any predicted radius of an ultracool dwarf. Discussion of this unphysical radius, or 
equivalently low luminosity is given in Sect.[5] 



tinct from field and young cluster L-dwarfs, with substantially 
different ratios of flux across the bands, a triangular shape to the 
H-band continuum, and a peak at longer wavelengths in the K- 
band. Spectral indices developed for cool objects give conflicting 
values, indicating the unusual nature of this object. The spectrum 
is best fit by a DUSTY model atmosphere with r ei f=1600K and 
log g- 4.5, though it is not yet possible to exclude the log g= 3.0 
atmosphere predicted by one formation model. We find that the 
apparent underluminosity of 2M1207 b is consistent with extinc- 
tion from large dust grains in a nearly edge-on disk, but more 
precise /-band and longer wavelength photometry and spec- 
troscopy will help to determine the nature of this important ob- 
ject. This spectrum of 2M1207b will serve as a key compari- 
son for upcoming near-infrared spectroscopy of planets imaged 
around stellar hosts. 
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